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ABSTRACT 

A geometrically nonlinear finite-element analysis of cohesive failure in 
typical joints is presented. Cracked-lap-shear joints were chosen for analysis. 
Results obtained from linear and nonlinear analysis show that nonlinear effects, 
due to large rotations, significantly affect the calculated mode I, crack 
opening, and mode II, inplane shear, strain-energy-release rates. The ratio of 
the mode I to mode II strain-energy-release rates was found to be 

strongly affected by the adhesive modulus and the adherend thickness. 
ratios between 0.2 and 0.8 can be obtained by varying adherend thickness and 
using either a single or double cracked-lap-shear specimen configuration. 

Debond growth rate data, together with the analysis, indicate that mode I strain- 
energy-release rate governs debond growth. Results from the present analysis 
agree well with experimentally measured joint opening displacements. 
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INTRODUCTION 


Joining structural components with adhesives offers many advantages over con- 
ventional mechanical fasteners: lower structural weight, lower fabrication cost, 

and improved damage tolerance have been demonstrated [1-5]. Fokker Aerospace has 
successfully built adhesively bonded aircraft for more than 30 years [4]. Signifi- 
cant new surface preparations and adhesives have appeared in the past 10 years that 
have enhanced the strength and durability of bonded joints. Even with 
all of the advantages and encouraging experiences with adhesive bonding, manufac- 
turers still hesitate to apply this technology to primary structural components, 
citing lack of confidence from poor understanding and concern for long term dura- 
bility. A Fokker report [5] showed that many of the past failures of bonded struc- 
tures involved poor design. Most often the poor design could be attributed to 
inadequate understanding of the adhesion failure mechanics. 

There are three dominant modes of mechanical failure in bonded joints: 

1) adherend failure (including delamination if one of the adherends is a laminated 
composite), 2) interfacial failure between adhesive and an adherend and 3) cohesive 
failure within the adhesive. The objective of the current study was to contribute 
to the understanding of the mechanics of cohesive failure. 

When a joint fails at the bond from cyclic loading, this failure mode is often 
called cyclic debonding-progressive separation of the adherends by cyclic failure 
of the adhesive. Previous researchers [6,7] have shown that the strain-energy- 
release rate defined from fracture mechanics principles can be used to correlate 
the cyclic debonding data. They used the total strain-energy-release rate as 
calculated from strength-of-materials principles. These results and more recent 
data [8] indicate that when mixed-mode loading is involved (i.e., peel and shear 
loads), cyclic debonding data can be correlated by an appropriate mix of mode I and 
mode II components of strain-energy-release rates. In general, these components 
of strain-energy cannot be calculated using strength of materials methods. A more 
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accurate analysis is needed so that the strain-energy-release components can be 
calculated from the local forces and displacements around the debond tip. This 

paper presents an analytical study of cohesive failure in adhesively bonded joints 
using a geometrically nonlinear two-dimensional finite element analysis. The analysis 
presented resolves the strain-energy-release rate at the debond tip into mode I and 
mode II components. 

Cracked-lap-shear (CLS) joints (Fig. la) were chosen for this study because 
they represent mixed-mode loading and large area bonds typical of many structural 
applications. They are also convenient specimens for laboratory tests on debond 
growth. The analytical study on these joints provides an insight into the geometric 
nonlinear effects. Furthermore, the analysis of these specimens illustrates the 
effect of adherend and adhesive material properties and joint dimensions on the 
mode I and mode II strain-energy-release rates. Data on debond growth in CLS 
specimens were analyzed to assess the relative influence of mode I and II strain- 
energy-release rates. The results of the present analysis are compared to those of 
other researchers. 
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NOMENCLATURE 

2 

area, m 

length of debond, m 

width of specimen, m 

material stiffness coefficients, GPa 

Young's modulus in and X2 directions, GPa 

Young's modulus in x-|,X2 directions for lap adherend, GPa 

Young's modulus in x.j,X2 directions for strap adherend, GPa 

Young's modulus of adhesive, GPa 

shear modulus, GPa 

shear modulus of adhesive, GPa 
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mode I strain-energy-release rate, J/m 

p 

mode II strain-energy-release rate, J/m 

total strain-energy-release rate (=>^ ^ 

element linear stiffness coefficients, kN/m 

element large displacement stiffness coefficients, kN/m 

element tangential stiffness coefficients, kN/m 

element geometric stiffness coefficients, kN/m 

length of lap adherend, m 

length of strap adherend, m 

applied load cycles 

applied load, kN 

applied load in direction at node A, kN 

thickness of adhesive, m 

thickness of lap adherend, m 

thickness of strap adherend, m 

strain energy, J 

displacement in x^. direction, m 

displacement in x^ direction at node 0, m 

debond opening displacement, m 

work done by external loads, J 

rectangular Cartesian coordinates, m 

strain components 

Poisson's ratio 

Poisson's ratio of adhesive 

total potential energy, J 

polynomial interpolation function for node 6 

stress components, GPa 
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residual force in direction at node A, kN 

Indices: 

i>j.m,n,r Cartesian directions, range=l,2 

A, e node numbers connecting an element, range=l,4 

BASIC JOINT CONFIGURATION 

Figure la shows the basic cracked-lap-shear specimen configuration chosen for 
the present analysis. This base specimen consists of a 203 mm long lap adherend 
bonded over a 254 mm long strap adherend. The strap adherend was 3.2 mm thick 
2024-T3 aluminum and the lap adherend was 1.6 mm thick unidirectional graphite/ 

epoxy. The adherends were bonded with a room temperature cure adhesive EA-934. 

The material properties of the adherends and adhesive were taken from the literature 
[9,10] and are shown in table 1. In the analytical results presented in this 
paper, one or more dimensions or properties were varied. Any variations from the 
base specimen are indicated wherever they occur. 

TABLE 1.- MATERIAL PROPERTIES. 


Material 

^1 

^2 

^21 

Gi2 


GPa 

GPa 

GPa 

Aluminum 

71.02 

71 .02 

0.33 

26.66 

Gr/Epoxy [9] 

130.93 

11.72 

0.0188 

4.46 

EA-934 [10] 

4.14 

4.14 

0.4 

1.475 


ANALYSIS 

The cracked-lap-shear joints were studied with a nonlinear two-dimensional 
finite element analysis. The analysis accounts for the kind of geometric non- 
linearity associated with large rotations in such unsymmetric joints. It was 
derived from an earlier analysis reported in reference 11. The computer program, in 

its present form, is called GAMNAS (Geometric ^nd Material Nonlinear Analysis of 
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structures). A two-dimensional idealization and the coordinate system used for the 
analysis are shown in fig. lb. The theoretical basis for the analysis is described 
in the following paragraphs. 


Governing Equations 

In the following development of governing equations Cartesian tensor notation 
is used. The governing nonlinear equations for individual elements can be derived 
using the principle of minimum total potential energy. The total potential energy, 
n, is given by equation (1) [12]. 

n--U-W = |/aijE,jdA-W (1) 

Geometric nonlinearity is included by using the nonlinear strain-displacement 
relation [13] in equation (1). 

ij 2 \ dx. 3x. 3x. 3x J 

The displacements u-j and within an element are approximated with interpolation 
functions, that is 

u. = 4>®u? (3) 

Using equations (1), (2), and (3) to calculate and minimize n with respect to u? 
yields 


Equations (4) are the set of governing nonlinear equations for the element. The 
equations for all elements are assembled to form a system of governing equations. 

The assembled nonlinear governing equations are solved iteratively using a 
modified Newton-Raphson method [12]. To implement this method, the assembled 
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tangential stiffness matrix is obtained by combining the individual element stiffness 
matrices. The element tangential stiffness matrix is obtained from equation (4) by 
partial differentiation as follows: 



Performing the differentiation and simplifying results in 
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^*^o^mn’ ^\^mn’ ^'^a^mn linear, large displacement, and geometric stiff- 

ness matrices, of each element, 'respectively. 

The elements used in the analysis were four-node, isoparametric quadrilaterals. 

To improve the performance of the elements in modeling bending deformation, a 
selective reduced numerical integration scheme [11,14,15] was used to evaluate the 

element stiffness matrices. The procedure for calculating the linear stiffness 
matrix (K^) is described in references 14 and 15. Application of reduced integration 
to nonlinear analysis is described in reference 11. 
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Calculation of Strain-Energy-Release Rates 
Strain-energy- release rates were calculated using a crack closure technique 
like that reported in reference 16. The forces transmitted through the node at the 
crack tip and the relative displacements of the two nodes on the crack boundary 
closest to the crack tip were used to calculate the energy required to close the 
crack. Because the rotations are significant, the nodal forces were calculated 
using equation (4). The forces calculated with equation (4) and the displacements 
were defined relative to a fixed coordinate system. These quantities must be 
transformed to obtain the components normal and tangential to the rotated bondline. 
Figure 2 shows the original and the transformed coordinates. The mode I strain- 
release-rate, , was calculated from the force and displacement components normal 
to the rotated bondline. n is due to in-plane shear mode of fracture and was 
obtained from the corresponding tangential components. The strain-energy-release 
rates and were obtained as (see fig. 2) 



(7) 


The total strain-energy-release rate,/^'(j. is the sum of and 

Finite Element Models 

The typical finite element model shown in figure lb has 880 nodes and 792 
elements. Other models had from 750 to 1100 nodes and 600 to 1000 elements and a 
convergence study showed that this level of refinement calculated the components 
of strain-energy-release rate within 3 percent of the values obtained with a much 
more refined mesh. Hence, this level of precision was deemed adequate. Boundary 
conditions are indicated in figure lb. Under load, because of the eccentric load 
path, the ends of the CLS specimen try to rotate. But the CIS specimen is normally 
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gripped such that the ends cannot rotate. So, to simulate grip loading, a multi- 
point constraint was imposed to prevent rotation of the loaded end of the model, 
i.e., all the u-j , displacements along the end are equal. For some of the configu- 
rations analyzed in this paper failure to properly model these grip conditions 
resulted in discrepancies up to 15 percent. 

RESULTS AND DISCUSSIONS 

In any structure, its response to applied load depends upon a number of 
geometric and material parameters. In this study, the cracked-lap-shear specimen 
was analyzed with the present nonlinear analysis, and the effect of several of 
these parameters on the mix of mode I and mode II components of strain-energy- 
release rate were calculated. First, the need for a geometrically nonlinear analysis 
is shown. Then the effects of varying the adherend dimensions and the implications 
of the dimensions on the CLS specimens are discussed. The effects of adhesive 
properties variations demonstrate the need to model the adhesive in the analysis. 

In past studies on adhesively bonded joints, simplified strength-of-material s 
analyses have been used to study strain-energy-release-rate components in CLS 
specimens. Results from the present analysis are compared with those from strength- 
of-materials approach. Finally, the present analytical results are compared with 
experimental results. 

The numerical results are presented for plane stress idealization of the 
joint. Plane strain idealization estimates lowered (up to 9 percent) values of 
mode I and mode II components of strain-energy-release rates, but the trends of 
results were the same as those for a plane stress idealization. 

Analytical Results 

Previously, researchers [17] showed that a geometrically nonlinear analysis 
was often needed when analyzing adhesively bonded structures. They showed that 
the stress distributions may differ significantly between a linear and nonlinear 
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analysis. In the current study of the CIS specimen the difference between the 
strain-energy-release rate components were calculated using geometrically nonlinear 
and linear analyses. The results differed and are shown in fig. 3 for a debond 
length of 50.8 mm. In this figure, is the strain-energy-release rate due to 
the opening mode of fracture, is the strain-energy-release rate due to the 

in-plane shear mode of fracture. The linear and nonlinear values of and 
differ noticeably with the differences widening as the load increases, as would be 
expected. The data indicate the desirability of a geometrically nonlinear analysis 
of CIS specimens. 

Analysis of the durability of bonded structures depends on the correlation of 
debond growth with the components of strain-energy-release rate. The basic CIS 
specimen configuration previously described was studied using the present nonlinear 
analysis. Figure 4 shows the mixed mode nature of the CLS specimen; the value of 
is approximately 50 percent of ^jj. For this particular configuration 
and are essentially constant over most of the lap length; this simplifies 
the correlation of strain-energy-release rates with debond growth data. 

In general, the strain-energy-release rates are not constant with debond 
length for specimens of all configurations. The adherend thicknesses relative to 
lap-length become important parameters in determining the variation of to 
ratio along the lap length. Figure 5 shows how vanes with debond length 

at two different values of lap adherend thickness. One of the cases, t^^ = 1.6 mm, 
corresponds to the base configuration as shown in fig. 4. The variation 
of with debond length is greatest for the specimen with thicker 

adherends. 

A further analysis was carried out to understand whether the range of the 
^l//^I ratio could be extended by changing the thickness of the lap adhered 
relative to the strap adherend. The results of this analysis, for a CLS specimen 
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with 50.8 nm debond, is shown in figure 6. Also shown is the trend for a configu- 
ration called a double-CLS specimen (DCLS). The DCLS specimen, which is symmetric 
as shown in figure 6, consists of an aluminum strap with two graphite/ epoxy lap 
adherends bonded to each side of the aluminum. Results show that a wide range of 
to ratios can be obtained by choosing the proper combination of adherend 
thicknesses and specimen configuration, i.e., CLS or DCLS. 

If each side of the DCLS debonds at the same rate, the specimen as a whole does 
not rotate. For this case the DCLS specimen has a higher to ratio than the 
CLS specimen (fig. 6). But previous tests [7] showed that symmetric debond growth 
rarely occurred. When the debond growth is unsymmetric, the ratio of to 
must be calculated for the unsymmetric configuration, and would not be as shown 
in fig. 6. 

In any structural analysis the validity of the results depends upon the 
accuracy of the material properties used. In the present analysis of the CLS 
specimen the influence of the shear modulus and Poisson's ratio of the adhesive on 
the strain-energy-release rate were studied. The results are shown in figure 7. 

For the configuration analyzed, the curves show that and noticeably 

affected by the adhesive shear modulus while Poisson's ratio does not appreciably 
affect either or /^jj. In fact, and are quite sensitive to changes in 
the adhesive shear modulus below 2 6Pa. This is significant because most current 
technology adhesives have shear moduli less than 2 GPa. Also, these results show 
the need to have accurate adhesive shear modulus data when calculating the strain- 
energy-release rates. 

Researchers [17,18] have shown that the through- the-thickness variation of the 
adhesive stresses is important. So, the influence of debond location on strain- 
energy- release-rates was studied. All the preceding analytical results in the 
report have been for a middle-of-bond cohesive failure. Calculations of strain- 
energy-release rates for a middle-of-bond cohesive failure and a failure mode that 
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is near the adhesive/aluminum interface are shown in figure 8. These results 
show that the calculated values of are greater for a failure that is near 
the interface, while those for are greater for a center of bond failure. 

Because of the expense of performing finite-element analysis, previous 
studies have often sought to simplify the analysis by considering a strength- 
of-materials approach [7,8]. Results derived from a strength-of-materials 
analysis for a CLS joint [8] are compared with finite-element results in 
figure 9. The strength of materials analysis was developed for metal adherends 
and did not account for the presence of adhesive. Therefore, comparative 
results were also obtained on a modified finite-element model for an aluminum- 
aluminum CLS joint which does not contain the adhesive. Figure 9 presents the 
versus the lap adherend thickness obtained with the present non- 
linear and strength-of-materials analyses. This comparison indicates that the 
two analyses agree only when the adherend thicknesses are nearly equal (t^^ = t^). 
Furthermore, considerable disagreement exists between the trends of the two 
analyses. These results suggest the strength of materials analysis is not appro- 
priate except for equal thickness adherends. 

Comparison of Analytical and Experimental Results 
Two examples of CLS joint tests are presented to show how well the present 
nonlinear analysis predicts test results and how it helps data interpretation. 

In the laboratory, only displacement and load can be measured directly. 

For the CLS specimen configuration (debond length = 50.8 mm) used in the study, 
the debond opening displacement, V^, was measured at the free end of the graphite/ 
epoxy laminate and is shown in figure 10. Excellent agreement was obtained 
between the present nonlinear analysis and the test data. The linear analysis 
overestimated the debond opening displacement, especially at the higher loads. 

Everett [19] hypothesized that is the principal mechanism in causing debond 
growth in a CLS specimen. Several tests were designed in reference 19 to show the 
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role of ^ during cyclic debond growth. For these tests it was hypothesized that 
if a sufficient clamping force was applied to a bonded joint ahead of the advancing 
debond, this clamping force would oppose the peel stresses that drive the debond and 
stop the debond at the clamp. 

A series of constant amplitude fatigue tests were run on the basic CLS specimen 
configuration under a maximum cyclic load per unit width of 482 kN/m. When a clamp- 
ing force per unit width of 9.6 kN /m or greater was applied to the specimen, the 
debond growth essentially stopped. The debond growth rate, da/dN, obtained from 
reference 19, is presented in figure 11. For the identical CLS specimen configura- 
tion and loading, the strain-energy- release rates were calculated. These results 
are shown in figure 11. Both and were nearly constant until the debond 
approached the clamp point (point D in figure 11). The peel stress component 
decreased to nearly zero close to this point. In contrast, increased near 

point D. These trends support the hypothesis that the peel stresses are likely 
to be the principal cause of cyclic debond growth, 

CONCLUSIONS 

A geometrically nonlinear finite element analysis was used to study the 
mechanics of cohesive failure in adhesively bonded structures. Cracked-lap-shear 
(CLS) specimens that experience the mixed-mode loading typical of structures were 
analyzed and tested. The results led to the following conclusions: 

1. Nonlinear effects due to large rotations are pronounced in CLS specimens 
and they significantly affect the calculated mode I and mode II strain-energy- 
release rates, and 

2. By varying adherend thickness and specimen configuration, the ratioy^j/ j 
can be varied widely. The length of the joint and the lap length significantly 
affect as the joint debonds. 

3. Adhesive properties significantly affect the calculated values of and 
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This indicates that the adhesive must be included in the analysis. 

4. The test results on debond growth imply that the mode I strain-energy- 
release rate is the main driver of cyclic debonding in the CIS specimen. 
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(a) Basic CLS specimen. 

Fig. 1 Basic configuration and finite-element model. 
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FINITE ELEMENT MESH 



(b) Finite-element mesh (X 2 coordinates are 10 times magnified). 
Fig. 1 Basic configuration and finite-element model. 



Fig. 2 Transformed coordinate system for C^colculotion. 
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Fig. 3 Comparison of and calculated from 

linear and nonlinear analysis. 
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Fig. 7 Effect of adhesive material properties on 
strain-energy-release rates. 
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Fig. 9 Comparison of results from strength-of -materials and 

present nonlinear analyses. 
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Fig. 10 Comparison of analytical results with experimental 
measurements of debond opening displacements v^j. 
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Fig. 11 Comparison of debond growth results with 
strain-energy-release rate components. 
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